Spin currents in a coherent exciton gas 
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Spin currents and spin textures are observed in a coherent gas of indirect excitons. Applied magnetic 
fields bend the spin current trajectories and transform patterns of linear polarization from helical 
to spiral and patterns of circular polarization from four-leaf to bell-like-with-inversion. 

Studies of electron spin currents in semiconductors led to the discoveries of the spin Hall effect [1- 6 , persistent spin 
helix [7], and spin drift, diffusion, and drag [SHTO]. There is also considerable interest in developing semiconductor 
(opto) electronic devices based on spin currents. An important role in spin current phenomena is played by spin-orbit 
(SO) coupling. It is the origin of the spin Hall effect and persistent spin helix. It also creates spin structures with 
the spin vector perpendicular to the momentum of the electrons in metals [11] and topological insulators p!2HT4] . 
While phenomena caused by SO coupling are ubiquitous in fermionic systems, they have yet to be explored in bosonic 
matter. Available experimental data for bosons include the optical spin Hall effect in photonic systems [TOl |l6l [T7] 
and spin patterns in atomic condensates [181 119]. Here, we report the observation of spin currents and associated rich 
variety of polarization patterns in a coherent gas of indirect excitons. Applied magnetic fields bend the spin current 
trajectories and transform patterns of linear polarization from helical to spiral and patterns of circular polarization 
from four-leaf to bell-like-with-inversion. We also present a theory of exciton transport with spin precession that 
reproduces the observed exciton polarization patterns and indicates trajectories of spin currents. 

Excitons - bound pairs of electrons and holes - form a model system to study spin currents of bosons [20j. SO 
coupling for an exciton originates from the combined Dresselhaus and Rashba effects for the electron and the hole 
[2TH2^ . An indirect exciton can be formed by an electron and a hole confined in separate quantum- well (QW) layers 
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FIG. 1: Polarization patterns in exciton emission, (a) Diagram of the CQW; e, electron; h, hole, (b) Schematic of 
exciton formation in the external ring (left) and LBS ring (right); Excitons (red) form on the boundary of hole-rich (blue) and 
electron-rich (green) areas. Exciton transport is indicated by red arrows, (c) A segment of the emission pattern of indirect 
excitons showing the external ring (left) and multiple LBS. (d,e) Patterns of linear Pun — (Ix — Iy)/{Ix + ly) (d) and circular 
Pa — — + I^-) (e) polarization of the emission of indirect excitons in the region shown in (c). Tbath = 0.1 K. 
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(Fig. la). The spatial separation reduces the overlap of electron and hole wavefunctions thus producing indirect 
excitons with lifetimes orders of magnitude longer than those of direct excitons [24l|25]. Due to their long lifetimes, 
indirect excitons can travel over large distances before recombination [26] and can cool down below the temperature 
of quantum degeneracy and form a condensate [TT. 

The condensation of indirect excitons was predicted to cause the suppression of exciton scattering [24] . The measured 
strong enhancement of the exciton coherence length [13] experimentally shows the suppression of exciton scattering. 
The suppression of scattering results in the suppression of the Dyakonov-Perel and Elliott- Yafet mechanisms of spin 
relaxation Furthermore, the spatial separation between an electron and a hole suppresses the Bir-Aronov-Pikus 
mechanism of spin relaxation for indirect excitons [29l [30] . The suppression of these mechanisms of spin relaxation 
results in a strong enhancement of the spin relaxation time in a condensate of indirect excitons. While the spin 
relaxation times of free electrons and holes can be short [29], the formation of a coherent gas of their bosonic pairs 
results in a strong enhancement of their spin relaxation times, facilitating long-range spin currents. 

Previous studies identified the external ring and localized bright spot (LBS) rings in the emission pattern of indirect 
excitons (Fig. Ic) as sources of cold excitons [13 . These rings form on the boundaries between electron-rich and 
hole-rich regions; the former is created by current through the structure (specifically, by the current filament at the 
LBS center in the case of the LBS ring), whereas the latter is created by optical excitation (Fig. lb), see ^3j and 
references therein. In [13 , we presented the studies of spontaneous coherence of indirect excitons in the ring region. 

Here, we present the studies of spin currents and associated spin patterns around the sources of cold excitons - the 
rings. We present the observation of patterns of circular polarization, corresponding to spin perpendicular to the QW 
plane. We also show that the observed polarization patterns are controlled by magnetic field: these data prove that 
the pattern of linear polarization corresponds to spin orientation rather than merely to the orientation of an exciton 
dipole. We also deduce trajectories of electron and hole spin currents from the measured polarization patterns. 

The exciton polarization currents and associated spin textures are revealed by the polarization pattern of the 
emitted light measured by polarization-resolved imaging (Fig. ld,e). Experiments are performed in an optical 
dilution refrigerator. The photoexcitation is nonresonant and spatially separated so that the exciton polarization is 
not induced by the pumping light. 
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FIG. 2: Control of polarization patterns. Measured (a) and simulated (b) patterns of linear polarization of the emission 
of indirect excitons Pun in the region of LBS for different magnetic fields perpendicular to the QW plane B. Measured (c) and 
simulated (d) azimuthal variation of Pun at a distance from the LBS center r = 9 /im for B = (black) and 7 T (red). Angles 
are measured from the ^-axis. Measured (e) and simulated (f) patterns of circular polarization of the emission of indirect 
excitons Pa in the region of LBS for different B. Measured (g) and simulated (h) cross sections of at x = for B = 
(black) and 7 T (red). For all data, Tbath = 0.1 K, the LBS is at (105,75) in Fig. Ic. See [31]. 
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The binding energy released at the exciton formation in the rings and the current filament at the LBS center heat 
the exciton gas. The former heating source depletes the exciton condensate in the rings [13^. The latter is so strong 
that no condensate forms in the LBS ring center and the exciton gas is classical there [13]. Excitons cool down with 
increasing distance r away from the heating sources so that they can approach the condensation temperature at r = ro 
where the condensation is detected by interferometric measurements [13^. 

The indirect excitons in GaAs CQW may have four spin projections on the z direction normal to the CQW plane: 
Jz = — 2,— l,+l,+2. The states = — 1 and +1 contribute to left- and right-circularly polarized emission and their 
coherent superposition to linearly polarized emission, whereas the states Jz = —2 and +2 are dark [29j. The electron 
and hole spin projections on the z-axis are given by J^, while in-plane projections of electron and hole spins can be 
deduced from the off-diagonal elements of the exciton spin density matrix, which can be obtained from the measured 
polarization pattern [3P. The exciton states linearly polarized along the axes of symmetry are generally split due to 
in-plane anisotropy induced by the crystallographic axis orientation and strain. 

The observed polarization patterns are qualitatively similar for both sources of cold excitons - the external ring 
and LBS ring. An LBS ring is close to a model radially symmetric source of excitons with a divergent (hedgehog) 
momentum distribution (Fig. lb) and we concentrate on the polarization textures around the LBS here. All LBS 
rings in the emission pattern show similar spin textures (Fig. ld,e). 

A ring of linear polarization is observed around each LBS center (Fig. Id, 2a). This ring is observed in the region 
r < ro where the exciton gas is classical. The linear polarization originates from the thermal distribution of excitons 
over the linearly polarized exciton states. Heating of the exciton gas by the current filament reduces the polarization 
degree in the LBS center and, as a result, leads to the appearance of a ring of linear polarization. No such polarization 
reduction is observed in the external ring, consistent with the absence of heating by current filaments in the external 
ring area (Fig. Id). 

A helical exciton polarization texture that winds by 27r around the origin emerges in the LBS area at r > ro 
where the condensate forms (the latter is measured by shift-interferometry ) , Fig. Id, 2a, 2c. The LBS exhibits a 
divergent hedgehog-shaped momentum distribution (Fig. 3d). The exciton polarization is perpendicular to exciton 
momentum (Fig. 2a, 3d). This produces vortices of linear polarization which emerge in concert with spontaneous 
coherence below the critical temperature [31j. The observed radial exciton polarization currents are associated with 
spin currents carried by electrons and holes bound into excitons as detailed below. 

Applied magnetic fields bend the spin current trajectories creating spiral patterns of linear polarization around the 
origin (Fig. 2a,c). The spiral direction of the exciton polarization current clearly deviates from the radial direction 
of the exciton density current (Fig. 2a,c). The control of the polarization patterns by magnetic field shows that they 
are associated with spin. 

Regular patterns are also observed in circular polarization (Fig. le, 2e). An LBS source of excitons generates a 
four- leaf pattern of circular polarization (Fig. le, 2e). This pattern vanishes with increasing temperature |31J. An 
applied magnetic field transforms the four-leaf pattern to a bell-like pattern of circular polarization with a strong 
circular polarization in the center and polarization inversion a few /im away from the center (Fig. 2e,g). 

Polarization patterns are also observed in the external ring region (Fig. ld,e). At low temperatures, the macroscop- 
ically ordered exciton state (MOES) forms in the external ring. The MOES is characterized by a spatially ordered 
array of higher-density beads and is a condensate in momentum space [13 . The polarization texture in the external 
ring region appears as the superposition of the polarization textures produced by the MOES beads with each being 
similar to the texture produced by an LBS (Fig. ld,e). A periodic array of beads in the MOES (Fig. Ic) creates 
periodic polarization textures (Fig. ld,e). The periodic polarization textures in the external ring region vanish above 
the critical temperature of the MOES [31 . 

Below we present a theoretical model which describes the appearance of the exciton polarization textures and links 
them to spin currents carried by electrons and holes bound into bright and dark exciton states. This model is based 
on ballistic exciton transport out of the LBS origin and coherent precession of spins of electrons and holes. The former 
originates from the suppression of scattering and the latter from the suppression of spin relaxation in the condensate 
of indirect excitons. The states with different spins are split due to the splitting of linearly polarized exciton states 
and SO interaction, which is described by the Dresselhaus Hamiltonian = f3e {k%(7x — kyCFy) for electrons and 
Hh = Ph {^x^x + ^y^y) holes [2TH23] (ke,h are electron and hole wave-vectors given by = kex^e/ {"^e + ^h)i 
kh = kex^h/i^e + ^/i), and rrih are in-plane effective masses of electron and heavy hole, respectively, kex is the 
exciton wave vector, Pe,h are constants, and (Jx^y are Pauli matrices). In the basis of four exciton states with spins 
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FIG. 3: Spin textures. Simulated in-plane exciton polarization (a), electron spin (b), and hole spin (c) patterns, (d) Schematic 
of exciton momentum (black arrows) and linear polarization (red lines) patterns. Schematic of effective magnetic fields given 
by the Dresselhaus SO interaction for electrons (e) and holes (f). Exciton polarization (g), electron spin (h), and hole spin (i) 
patterns in applied magnetic field B — 7 T. The lines (arrows) and the color visualize the orientation of the linear polarization 
(spin). 



Jz = +1,-1, +2, —2, the coherent spin dynamics in the system is governed by a model matrix Hamiltonian: 
■ - {gn - 9e)^JiBB/2 -8^ kepee''^ khf^ne-'^ 

kei^ee'^ khPhe-'"^ Ed - {gn + 9e)^BB/2 -Sd • U 

khi^he'^ kePee-'^ Sd Ed + {gn + 9e)^BB/2 _ 

where £^5 and Ed are energies of bright and dark excitons in an ideal isotropic QW, and 5d describe the effect 
of in-plane anisotropy resulting in the splitting of exciton states linearly polarized along the axes of symmetry. The 
angle if is measured from the x axis. The details of this model are presented in [31 . Exciton propagation out of the 
origin governed by this Hamiltonian results in the appearance of a vortex of linear polarization with the polarization 
perpendicular to the radial direction and a four-leaf pattern of circular polarization in 5 = 0, as well as spiral patterns 
of linear polarization and bell-like patterns of circular polarization in finite magnetic fields. This model qualitatively 
reproduces the main features of the experiment for both linear (Fig. 2a-d) and circular (Fig. 2e-h) polarizations. 

This model describes the exciton polarization currents and allows deducing the spin currents carried by electrons 
and holes bound to excitons as detailed in [31 . Figure 3b, c show the electron and hole spin textures deduced from the 
measured exciton polarization texture (Fig. 2a, 3a). One can see that both the electron and hole spin tend to align 
along the effective magnetic fields given by the Dresselhaus SO interaction Bejj(e) = -^^jf^{—ke,x-,ke,y)-, ^eff{h) = 
Qhi^B (^h,x^ kh^y) (Fig. 3e,f), consistent with the model. The patterns of Pun corresponding to the simulations in Fig. 
3 are shown in Fig. 2b. 

The model can be improved by including non-linear effects. In [31], we present simulations of exciton spin currents 
using Gross-Pitaevskii type equations, which treat the excitons as a coherent field outside the LBS center and include 
dispersion and interaction. The simulation results are similar to that within the density matrix approach and are in 
agreement with the experiment. Non-linear spin-related phenomena form interesting perspectives for future studies. 
In conclusion, long-range spin currents governed by spin-orbit interaction and controlled by applied magnetic field 
have been observed in a coherent exciton gas. 
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THEORY OF EXCITON SPIN CURRENTS: DENSITY MATRIX APPROACH 
Ballistic exciton transport with coherent spin precession 



In zinc-blend semiconductor quantum wells (e.g. in GaAs/AlGaAs system), the lowest energy exciton states are 
formed by electrons with spin projections to the structure axis of +1/2 and —1/2 and heavy holes whose quasi- spin 
(sum of spin and orbital momentum) projection to the structure axis is +3/2 or —3/2. Consequently, the exciton 
spin defined as a sum of electron spin and heavy hole quasi-spin may have one of four projections to the structure 
axis: — 2,— l,+l,+2 [1 . These states are usually nearly degenerate, while there may be some splitting between 
them due to the short and long-range exchange interactions. Here we derive the exciton Hamiltonian in the basis 
of +1, — 1,+2, —2 states, accounting for the spin-orbit interaction [2|l3l, long- and short-range exchange interactions 
[4] and Zeeman effect, but neglecting exciton-exciton interactions. We consider excitons propagating ballistically in 
plane of a quantum well. We shall characterize them by a fixed wave- vector kex- 

In order to build the 4x4 matrix Hamiltonian for excitons, we start with simpler 2x2 Hamiltonians describing 
the spin-orbit Dresselhaus effect and Zeeman effect for electrons and holes. 

The electron Hamiltonian in the basis of (+1/2, —1/2) spin states writes: 

He = f^e{ke,xCrx - ^e^yCTy) - ]^gel^BB(7z. (1) 

Here is the electron ^-factor, /i^ is the Bohr magneton, 5 is a magnetic field normal to the quantum well plane, 
Pe is the Dresselhaus constant describing spin-orbit interactions of electrons, the Pauli matrices are 

Hence 





'0 1" 




'0 -z" 




■ 1 




1 


^ = 


i 




-1 



— \gel^BB 



e,y 



/^e(^e,x + ^^e,y^ 



i(f) 



(2) 



Here (j> is the angle between the electron wave-vector kg and the chosen x-axis. The exciton Hamiltonian needs to be 
written in the basis of (+1, —1, +2, —2) exciton states, which correspond to (—1/2, +1/2, +1/2, —1/2) electron states. 
The electron spin-flip couples +1 and +2 states and —1 and —2 states. For each of these two couples of states we 
apply the Hamiltonian (2), which results in the following electronic contribution to the 4x4 exciton Hamiltonian: 



ge^iBB/2 

-9el^BB/2 

kePee'* 

fce/3ee-'* 



keAe-'* 

kep.e"^ 

-9eliBB/2 





(3) 



The similar reasoning applies to the heavy hole contribution to the Hamiltonian: The hole Hamiltonian written in 
the basis of (+3/2, —3/2) states is 

Hh = Ph{kh,xcrx + kh,ycry) - ]^ghl^BB(7z. 
Here g^ is the heavy hole ^-factor, {3^ is the Dresselhaus constant for heavy holes [21 [3]. Hence, 



-\ghl~iBB (3h{kh,x - ^kh^y) 



Ph{kh,x + ^^/i,2y) 



Phkhc'^ 



^gh/J^sB 



(4) 



^ghj^sB 

Here is the angle between the hole wave vector kh and the chosen x-axis. The exciton Hamiltonian is written in the 
basis of (+1, —1, +2, —2) exciton states, which correspond to (+3/2, —3/2, +3/2, —3/2) hole states. The hole spin-fiip 
couples +1 and —2 states and —1 and +2 states. For each of these two couples of states we apply the Hamiltonian 
(4), which results in the following hole contribution to the 4x4 exciton Hamiltonian: 



-9h^iBB/2 

ghl^BB/2 

khPhS-''^ 
khPhe'^ 



i(f) 



khPhe 
khPhe''^ 
-ghl^BB/2 

gh^BB/2 



(5) 
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For the translational motion of an exciton as a whole particle the exciton momentum is given by 
Vex = (^e + "^hhj^exi where TTig and rrihh are in-plane effective masses of an electron and of a heavy hole, 
respectively, Vex is the exciton speed. Having in mind that the exciton translational momentum is a sum of electron 
and hole translational momenta given by Pe,/i = ^e,/i/iVe,/i, Ve,/i being the electron (hole) speed, one can easily see 



that 



Having in mind that Pg 



hkp h we have kp 



me+ruhh 



me+mhh 



kp 



Besides the contributions from electron and hole spin orbit interactions and Zeeman splitting, there may be a purely 
excitonic contribution to the Hamiltonian, which is composed from the Hamiltonian for bright excitons written in the 
basis (+1, —1): 



Hh — E})I — S^ax 



-Sb Eb 



and the Hamiltonian for dark excitons written in the basis (+2, —2): 

El^d = EdI 



Ed -Sd 
-Sd Ed 



(6) 



(7) 



where / is the identity matrix. The terms with and Sd describe the splittings of bright and dark states polarized 
along X and y axes in the plane of the structure due to the long-range exchange interaction. Such splitting may 
appear due to some in-plane anisotropy in the structure induced by strain or monolayer fluctuations of interfaces. We 
assume that it induced the splitting of X- and F-polarized excitons, while it can be easily generalized to the splitting 
in diagonal or random axes. E^ — Ed is the splitting between bright +1 and —1 and dark +2 and —2 exciton states 
which may be also split due to the short range exchange interaction. Note that Eqs. (6,7) can be simply obtained 
from the exciton Hamiltonians written in the basis of (X, Y) polarizations. E.g. for the bright excitons: 



HxY 



Eb-h 

^6 + ^6 



(8) 



where C 



1 

V2 



1 1 

i —i 



1 

V2 



C ^ HxyC 



are the transformation matrices from linear to circular basis and vice 



versa [5]. The same reasoning applies to the dark excitons as well. The sum of Hamiltonians (6,7) in the 4x4 basis 
writes: 



Ho 





-5b 








-Sb 


Eb 














Ed 


-Sd 








-Sd 


Ed 



(9) 



Let us consider the excitons propagating with a wavevector kg^^. We shall describe them by a spin density matrix 
p = 1^ >< ^1, where ^ = (^+i, ^+2, ^-2) is the exciton wave-function projected to four spin states. The 
elements of this density matrix pij are dependent on the distance from the excitation spot r = Vpxt and the polar 
angle 0. The elements of the upper left quarter of the density matrix are linked to the intensity of light emitted by 
bright exciton states I and to the components of the Stokes vector Sx, Sy and Sz of the emitted light: 



Pii — 2 + /^i2 — Sx — iSy^ P21 — Sx ^ iSy^ P22 — 2 



These expressions can be summarized using the Pauli matrices as 



Pii P12 
P21 P22 



Sz. 



(10) 



S(j, where / is the identity 



matrix. 

The components of the Stokes vector are directly proportional to the polarization degree of light measured in XY 
axes, diagonal axes and the circular basis. The circular polarization degree of light emitted by propagating excitons 
can be obtained as 



Pc = 2Sz = {pU - p22)/{pll +P22), 



(11) 
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the linear polarization degree can be found from 

pi = 2Sx = {pi2 + P21)/ {pii + P22), 



(12) 



the linear polarization degree measured in the diagonal axes (also referred to as a diagonal polarization degree) is 
given by 



Pd = '^Sy = i{pi2 - P21)/ {Pii + P22)' 
The dynamics of this density matrix is given by the quantum Liouville equation: 

^n'i = [H,p], 



(13) 



(14) 



where the Hamiltonian is composed from the electron, hole and exciton contributions given by Eqs. (3,5,9) as follows: 



H = 



{9h - 9e)l^BB/2 

-Sb Eb 



{9h - ge)pBB/2 
khPhe-'^ Ed 
kePee-'* 



fee/See"** 

khPhe''^ 
{Oh + 5e)MB-B/2 
-5d Ed 



khPhS-''^ 

kePee'* 

-Sd 

{gh + ge)i^BB/2 



(15) 



The Hamiltonian (15) includes the electron, hole, and exciton contributions. Magnetic field affects the electron and 
hole contributions via the Zeeman splitting. Its effect on the exciton contributions, originating from the change of 
short and long-range exchange interactions, is neglected in the model. A qualitative agreement with the experiment 
justifies this approximation. 

Note, that in a similar way one can describe the Rashba effect for electrons and holes on the exciton spin density 
matrix. Our estimations show that for the value of bias we use in these experiments the Rashba effect is much 
weaker than the Dresselhaus effect [6 . Therefore we limit ourself to the consideration of the Dresselhaus effect for 
electrons and holes. In order to make sure that the observed exciton polarization textures are indeed governed by 
the Dresselhaus effect, we have performed also the simulations accounting for the Rashba instead of Dresselhaus 
mechanism of spin-orbit coupling. These simulations produce exciton polarization patterns qualitatively different 
from the experimental data. 



An example of how the Dresselhaus effect affects the polarization of propagating exciton 

In order to obtain the spatial distribution of Stokes vector components in the cw regime we assume that all excitons 
propagate in radial directions from a point-like or a ring-like source. Their polarization state in a point characterized 
by the polar coordinates (r, (j)) is readily obtained from the elements of the density matrix p(t, 0) with t = r/vex- The 
exciton speed Vex governs the spatial scale of the polarization textures. 

Let us consider the simplest example of how the Dresselhaus effect affects the polarization of propagating excitons. 
In order to do it, we shall commute both parts in Eq. (14) with the Hamiltonian. As a result we shall have: 



^n'^ = mH,p]]. 

Now we take a time derivative from both parts of Eq. (14) and substitute the expression (16) in its right part: 



(16) 



h'^ = [H,[H,p]]. 



(17) 



Let us suppose that initially we have an exciton state composed by bright excitons linearly polarized along x-axis and 
dark excitons linearly polarized along y-ax.is. It is described by the density matrix: 



Po 



110 
110 
1-1 
0-11 
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How this matrix would evolve in time due to the Dresselhaus effect on the electron spin? Let us assume B = for 
simplicity, and calculate the commutator of the Hamiltonian (3) 
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— cos 










— cos 
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The double commutator in the right part of Eq. (14) can be now found as: 

[Fe,[Fe,P0]]=4(Me)' 



COS 0e 

cos(/)e^^ 











COS (j)e~ 






— cos 0e*' 




Let us substitute this expression back to Eq. (17) and look at the dynamics of the element pi2 = Sx — iSy 
describing the linear and diagonal polarization degrees of bright excitons. One can easily see that 
^^^2 ^ = —fl^ ^ ^^dt^^^^^ ~ ^{kePe)'^ cos(l)e~'^^ . Separating the real and imaginary parts of this equation, we 
find the dynamics of the Stokes vector components: h^^-^^ = -4:{kePef ^os'^(j)\h^^-^^ = -2(A:e/5e)^ sin20. 



Our initial conditions are: dit t = 0, Sx = Sxo > 0, Sy = 0. One can see that ior (p = — 7r/2,7r/2 the polarization does 
not change: Sx = Sxo^Sy = 0. For ^ = 0,7r, we have 5*^ = 0, while Sx decreases and, eventually, inverts its sign. The 
negative Sx corresponds to F-polarization. For any other value of both Sx and Sy change with time. Namely, Sx 
decreases, Sy builds up. The sign of Sy is negative (corresponds to the polarization along (1,-1) axis) ifO<0<7r/2 
and TT < (j) < 37r/2. The sign of 5*^ is positive (corresponds to the polarization along (1,1) axis) if7r/2<0<7r and 
37r/2 < < 27r. This describes a polarization vortex. 

In the same way one can show that the Hamiltonian (9) describing the exchange induced linear polarization splitting 
converts any linear polarization different from X and Y polarizations to the circular polarization. 



Spin currents carried by electrons and holes bound into excitons 



It is important to note that the present formalism addresses the spin part of the exciton wavefunction, which is a 
product of electron and hole spin functions. E.g. the probability to find the exciton in the spin state +1 is given by 
a product of probabilities to find an electron in the spin state —1/2 and the heavy hole in the spin state +3/2. The 
four component exciton wave-function 



^ = (^+1, ^+2, ^-2, ) = (^e,-l/2^/.,+3/2, ^e,+l/2^/^,-3/2 , ^e, + l/2^/^,+3/2 , ^e,-l/2^/.,-3/2) , 



(18) 



where ^e,+i/2 and ^e,-i/2 are the components of the electron spinor wave function, '^h,-\-3/2 and '^h,-3/2 are the 
components of the heavy hole spinor wave function. We shall normalize exciton, electron and hole wave functions to 
1, namely: 



^+1"^+! + ^-l^li + ^+2^+2 + ^-2^12 = ^e,+ |^e,+ | + "^e-^K-^ 

Now, the exciton spin density matrix is given by 



h,- 



2 ^^5— 2 



(19) 



4'+i*;2 *+i*-2 

*_i*+2 *-l*-2 

*+2*;i *+2*-l *+2*;2 *+2*-2 

*_2*+i *-2*^i *-2*+2 *-2*-2 
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e 



in 



1 ^* 1 ^, • 

1 ^* 1 ^, ' 



e,+ 



" e,+ 



(20) 

This representation allows us to obtain useful links between the elements of exciton, electron and hole density matrices, 
in particular: 



e,+ 



1^*1 ^ ^i^*_i 



P22 + P33 P24 + P31 



(21) 



Pll + P33 Pl4 + P32 
P23 + P41 P22 + P44 



(22) 



We know that the components of electron and hole density matrices are linked with the projections of electron and 
hole spins as: 



2 



Sh,z Sh,x '^Sh,y 



Sh,z 



(23) 



(here, for the heavy hole we have assigned spin +1/2 to the state +3/2 and spin —1/2 to the state —3/2 accounting 
for the orbital momentum of these states of +1 and —1, respectively). Therefore, the z-component of spin polarization 
carried by electrons can be expressed as 



Se,z = {P22 + P33 - Pll - P44)/2, 

and the 2;-component of spin polarization carried by holes can be expressed as 

Sh,z = {Pll + P33 - P22 - P44)/2. 



(24) 



(25) 



The in-plane component of electron and hole spins can be extracted from the off-diagonal elements of the density 
matrix. Namely, the x-component of electron spin is given by 



Se,x = {Pl3 + P31 + P24 + P42)/2, 

the x-component of the hole spin is given by 

Sh,x = (Pl4 + P23 + P32 + P4l)/2, 

the ^/-component of electron spin is given by 

Se,y = i{-Pl3 + P31 + ^24 " P42)/2, 

and the ^/-component of the hole spin is given by 

Sh,y = KP14. - P23 + P32 " P4l)/2. 



(26) 
(27) 
(28) 
(29) 



Effective fields 



Spin-orbit effect for electrons and holes can be accounted by introducing an in-plane effective magnetic field: 

1 



effich) = -^ge,hPB (Be//(3-) . 



In the case of Dresselhaus Hamiltonian for electrons: 

1 



gePB'Beff = Pe (ke,x, -^e,y) • 



For holes the field is collinear to the wave- vector: 



- ^dhPB'^eff = Ph {kh,x, kh,y) . 



(30) 
(31) 
(32) 
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FIG. SI: Simulated polarization patterns. Simulated in-plane exciton polarization in magnetic field B = (a) and 7 T 
(b). The lines and the color visualize the orientation of the linear polarization. 



Simulations 



The system is modeled using Eq. (15) and parameters are adjusted to match the experimental data (see Fig. 2 in 
the main text). The same set of parameters is used for linear and circular polarizations and for all magnetic fields. 
The parameters are /3e = 2.7 /ieV/im, /^/^ = 0.92 /ieV/im, ^5 = 0.5 /ieV, (5^^ = — 13 /ieV, Ei^ — = 5 /ieV, kex = 15.4 
/im~^, T = 0.1 K, = —0.01, and gh = +0.0085. Some of the parameters were obtained in earlier studies and some 
of them were put as fitting parameters. All fitting parameters were chosen to be consistent with the data published 
in earlier studies: Pe is consistent with our earlier measurements [6], - with the calculated values |2l|^, T - with 
the calculated temperature of indirect excitons for the studied structure [T, exciton splittings - with typical exciton 
splittings in GaAs structures [8] , ^f-factors - with typical ^-factors in quantum wells with ~ 8 nm width [9 , and kex 
was taken within the light cone in the structure (we checked that the model leads to qualitatively similar patterns for 
various kex)- 

The images are presented with spatial averaging over 1.5 /im. The initial exciton state considered by this model 
is a ring around the LBS center where the exciton gas is classical. The ring radius is taken 4 jam. On this ring, the 
simulations consider the classical exciton energy distribution with T = 0.1 K, kex = 0- No simulations were performed 
inside this ring, the polarization in the ring is shown there. Beyond this ring, the exciton gas is coherent and the 
simulations consider ballistic exciton transport with coherent spin precession. 

Figure SI (an enlarged version of Fig. 3a, g in the main text) presents simulated patterns of in-plane projection 
of the Stokes vector of emitted light which directly maps the pseudospin of bright excitons (see e.g. [101 [H]. The 
patterns of light polarization corresponding to these simulations are shown in Fig. 2b in the main text. Figure SI 
visualizes a vortex (a) and a spiral (b) pattern of linear polarization in zero (a) and a finite (b) magnetic field. The 
simulations produce the observed exciton polarization textures as described in the main text. 



THEORY OF EXCITON SPIN CURRENTS: GROSS-PITAEVSKII EQUATIONS 

The spin density matrix theory is convenient for the description of (partially) coherent and (partially) polarized 
exciton gases. However, the treatment of non-linear effects in a partially coherent system is a challenging task. An 
approach, which we consider in this section is to assume a perfectly coherent condensate of excitons. In this case, the 
excitons can be described by a spatially dependent four-component wavefunction, (?/;+i(x), ?/^_i(x), ?/^+2(x), V^_2(x)). 
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a B=OT B=1T B=3T B=5T B=7T 




X (yc/m) 

FIG. S2: Gross-Pit aevskii simulation. Spatial distribution of the linear (a) and circular (b) polarization degrees vs magnetic 
field. As in the case of the density matrix calculations, the results are spatially averaged over 1.5/im. Parameters: ruex = 0.21me, 
/3e = lO/xeV/im; I3h = /^e/3; St = 2/ieV; 6d = 0.5/ieV; Eh — Ed = 2/ieV; W = 0.2a, an = 1/ieV. The source area was taken 
circular with a radius of l//m. 



The dynamics of the wavefunction can be described by the Gross-Pitaevskii equation: 



ih- 



dt 



/V^+i(x)\ 
V^-i(x) 

V^+2(x) 

vv^-2(x)y 



H 



( ^+i(x) \ 




/^+i(x)\ 




/^*l(x)tA+2(x)^- 


-2(X) \ 


^/'-i(x) 


+ an(x) 


^/'-i(x) 


+ 


V';i(x)V'+2(x)^/'- 


-2(X) 


^/'+2(x) 


V'+2(x) 


V'-2(x)V'+l(x)^/'- 


-i(x) 


V ^/'-2(X) y 




\^-2(x) y 




V V'+2(x)V'+l(x)^/'- 


-i(x)/ 



(33) 



Here H is the same Hamiltonian (15), although to correctly define the operation of the /c— dependent spin-orbit 
terms one should replace ke = — — kex and kh = — ^ — kex where kex = — ^V. Note that the Gross-Pitaevskii 

^ ^ me-\-mhh me-\-mhh 

equation allows to work with a distribution of wavevectors, and accounts for the dispersion of excitons via the term 
^ ^ where rriex is the exciton effective mass. The last two terms in Eq. 



2me 



33 



are nonlinear terms. In the initial 

simplified approach presented here we assume a spin-independent scattering rate^a, where each spin fraction scatters 
with the total density, n(x) = |?/^+i(x)p + ?/^_i(x)^ + |?/^+2(x)p + |?/^_2(x)p. The last term, proportional to W, 
represents a parametric scattering process unique to indirect exciton systems where two bright excitons convert into 
two dark excitons (or vice versa). Note that all non-linear scattering processes conserve the total spin projection of 
excitons. 

To describe the excitation of our system, we introduce a fixed wavefunction boundary condition along a circular 
boundary representing the edges of the classical region around the LBS center. We choose a linearly polarized dark 
exciton density at this boundary, assuming that dark excitons have lower energy than bright excitons. We do not aim 
to use the Gross-Pitaevskii equation in the hot LBS center where the exiton gas is classical, however, we expect it to 
offer a qualitative description of the propagation of coherent excitons away from the LBS source in the region beyond 
the LBS center where the exciton gas is coherent. The indirect excitons have a long lifetime and thus we employ an 
absorbing boundary condition to allow solution of the problem. The absorbing boundary introduces a loss mechanism 
in the system, such that a steady state is reached where the excitons excited at the LBS balance the flow of excitons 
away from the studied region (we assume that upon crossing the absorbing boundary excitons never return). 

By solving Eq. [33] numerically for the steady state, the linear and circular polarization degrees are calculated for 
increasing magnetic fields (Fig. S2). The simulated patterns are qualitatively similar to both the experimentally 
measured patterns and the patterns simulated using the spin density matrix theory. 
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EXPERIMENT 
Experimental setup 



optical dilution refrigerator 



sample 



A/4 A/2 



detector 



linear 
polarizer 



FIG. S3: Experimental setup. Schematic of polarization-resolved imaging. 



The schematic of the polarization-resolved imaging experiment is presented in Fig. S3. The sample is in an optical 
dilution refrigerator. Light emitted by the sample is collected and made parallel by an objective inside the refrigerator. 
Polarization selection is done by a combination of quarter- wave plate (A/4), half- wave plate (A/2), and linear polarizer. 
The linear polarizer is aligned such that y-axis polarized emission is transmitted to the detector. The detector is a 
combination of an interference filter of linewidth ±5 nm adjusted to the emission wavelength of indirect excitons 
A = 800 nm (due to the interference filter, only the emission of indirect excitons is measured: the contribution of the 
weak emission of direct excitons or any other emission, such as low-energy bulk emission, is cut off by the interference 
filter), a spectrometer operating in dispersionless zeroth-order mode, and a liquid nitrogen cooled CCD. 

Measurements of linear polarization. The quarter- wave plate is aligned so that the fast and slow axis are along the 
X- and y-dix.is. Hence, x- and ^/-polarized emission and ly are transmitted unchanged. To measure the half- wave 
plate fast axis is aligned along the ?/-axis. Then, ly is transmitted unchanged, and is then transmitted through the 
linear polarizer. To measure Ix, the half- wave plate axis is aligned 45° relative to the y-axis. Then, Ix is rotated to 
the 7/-axis and is transmitted through the linear polarizer. 

Measurements of circular polarizations. The quarter- wave plate is aligned so that the fast and slow axis are rotated 
by 45° with respect to the ?/-axis. Then, circularly polarized emission /cr+ and 1^- are converted to x- and ^/-polarized 
light. This light is then selected as described above. 

In all experiments, the photoexcitation is nonresonant (> 400 meV above the energy of indirect excitons) and 
spatially separated (the 10 /im-wide excitation spot is > 80 /im away from both the LBS and external ring) so that 
neither the exciton polarization nor coherence is induced by the pumping light. 



Patterns of linear polarization and coherence in the LBS ring region 

The ring of linear polarization vanishes with increasing temperature (Fig. S4d-f). The vortex of linear polarization 
vanishes with increasing temperature (Fig. S4d-f,j). 

The coherence of an exciton gas is imprinted on the coherence of emission, which is described by the first-order 
coherence function gi{Sx). In turn, this function is given by the amplitude of the interference fringes Ainterf(^x) in 
'the ideal experiment' with perfect spatial resolution. In real experiments, the measured Ainterf(^^x) is given by the 
convolution of gi{Sx) with the point-spread function (PSF) of the optical system used in the experiment p2|- Both 
for a classical gas and quantum gas gi{Sx) is close to 1 ai Sx = and drops with increasing Sx within the coherence 
length ^. The difference between the classical and quantum gas is in the value of ^. For a classical gas, (_ci is close 

to the thermal de Broglie wavelength A^^^ = y^^p^, which is well below the PSF width in the studied temperature 
range (^d@o.iK ^ 0-3 /im, the PSF width is ~ 1.5 /im). Spontaneous exciton coherence with a large coherence length 
^ is observed in the region of the polarization vortex (Fig. S4g, black squares). In contrast, ^ is short in the region 
of the polarization ring (Fig. S4g, red circles), which is close to the hot LBS center. Large ^ ^ in the region of 
the polarization vortex indicates a coherent exciton state with a much narrower than classical exciton distribution in 
momentum space, characteristic of a condensate. Small ^ ~ 1.5 jam in the region of the polarization ring indicates a 
classical exciton state with ^ measuring the PSF width. The coherence measurements are discussed in detail in ^T3] . 
With reducing temperature, the exciton spin textures emerge in concert with coherence (Fig. S4h and S4i). 
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FIG. S4: Pattern of linear polarization in the LBS region: temperature dependence. Measured PL intensity in the 
region of LBS centered at (105, 75) in Fig. Ic at Tbath = 0.1 (a) and 7 (b) K. (c) PL Intensity profile through the center of the LBS 
for Tbath — 0.1 (black) and 7 (red) K. Measured linear polarization of the emission of indirect excitons Pun — (Ix — ly) / (Ix -\- ly) 
at Tbath =0.1 (d) and 7 (e) K. Pun profile (f) vs x through the LBS center and (j) vs azimuthal angle measured from the 
y-axis at r = 9 /xm from the LBS center for Tbath = 0.1 (black) and 7 (red) K. (G) The exciton coherence degree measured by 
shift-interferometry: Interference visibility Ainterf vs. shift Sx for the vortex of linear polarization (squares), 18 /xm left of LBS 
center, and the polarization ring (points), 2 /am left of LBS center for the LBS centered at (80, 105) in Fig. Ic. Tbath = 0.1 K. 
Emergence of the vortex of linear polarization and spontaneous coherence at low temperatures: (h) The amplitude of azimuthal 
variation of Pun at r 10 fim from the LBS center and (i) interference visibility Ainterf in the polarization vortex for Sx — 2 
/xm vs. temperature. 



Patterns of circular polarization in the LBS ring region 



The four-leaf pattern of circular polarization vanishes with increasing temperature (Fig. S5a-c,e). At low tempera- 
ture, the flux of excitons from the LBS center is anisotropic: the emission intensity is enhanced along the polarization 
direction in the polarization ring (Fig. S6a). To quantify the flux anisotropy, we use the ratio of maximum to min- 
imum in the azimuthal variation of the total emission intensity of indirect excitons Imax/Imin at distance r = 8 /im 
from the origin. At high temperature, the exciton flux anisotropy vanishes (Fig. S5d). The four- leaf pattern and flux 
anisotropy emerge in concert (Fig. S5d and S5e). The four- leaf pattern of circular polarization is associated with a 
skew of the exciton fluxes in orthogonal circular polarizations (Fig. S6). 
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FIG. S5: Pattern of circular polarization in the LBS region: temperature dependence. Measured circular polar- 
ization of the emission of indirect excitons Pa = — + /^-) in the LBS region at Tbath = 0.1 (a) and 7 (b) K. 
(c) The azimuthal variation of P^r at r = 9 /xm from the LBS center for Tbath = 0.1 (black) and 7 (red) K. (d) The ratio of 
maximum to minimum in the azimuthal variation of the total emission intensity of indirect excitons Imax/ Imin at r = 8/im 
from the LBS center vs. temperature, (e) The amplitude of variation of Pa around the LBS centered at (105,75) in Fig. Ic, 
vs. temperature. 




FIG. S6: Pattern of circular polarization in the LBS region: azimuthal dependence, (a) Azimuthal variation of the 
emission intensity of indirect excitons at r = 8 /xm in (blue) and a~ (black) polarizations. Angles are measured from the 
y-axis. (b,c) Traces of the (red) and a~ (green) emission peak around = 270° [see (a)]; The emission image in cr^ (b) and 
(T~ (c) polarization is also shown. 



Polarization patterns in the external ring region 



The fragmentation of the external ring vanishes with increasing temperature (Fig. S7a,d,g), consistent with earlier 
results [14]. The periodic polarization textures in the region of the external ring vanish with increasing temperature 
(Fig. S7b,c,e,f,h,i,k,l). 

Both the position of the external ring and the wavelength of the exciton density wave formed in the ring are 
controlled by the laser excitation indicating that the exciton density modulation in the MOES is not governed by 
defects in the sample. The observation of polarization textures around the MOES beads, which are not associated 
with defects, shows that the polarization textures do not arise due to defects. 



16 




FIG. S7: Polarization textures in the region of external ring: temperature dependence. Polarization textures in 
the region of external ring. (a,d) Emission of indirect excitons in the external ring region. Pattern of linear (b,e) and circular 
(c,f) polarization in the region of external ring. Tbath = 0.1 (a-c) and 7 (d-f) K. PL intensity (g), Pun (h), and Pa (i) profiles 
taken along the external ring (along an axis 9 degrees of vertical positioned at maximum variation) for Tbath = 0.1 (black) and 
7 (red) K. PL intensity (j), Pun (k), and P^ (1) profiles taken across the external ring (along the x-axis positioned at maximum 
variation) for Tbath = 0.1 (black) and 7 (red) K. The image in (j) shows the external ring at Tbath = 7 K. 



The role of excitonic effects in the formation of spin textures in the exciton system 

The excitonic terms describing the spin precession in the Hamiltonian due to splitting of exciton states are present 
for excitons and absent for free electrons and holes. To verify the role of excitonic effects in the formation of spin 
textures, we performed the simulations with and without the excitonic terms in the Hamiltonian (15). No spin 
texture forms in the absence of excitonic terms, see Fig. S8. This shows that the excitonic effects are important for 
the formation of spin textures in the exciton system. 
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FIG. S8: Importance of the excitonic effects. Simulated Pun and Pa (a,b) with the excitonic effects [nonzero excitonic 
terms in Hamiltonian (15), the values are presented in Section Simulations] and (c,d) without excitonic effects [zero excitonic 
terms in Hamiltonian (15), 6b = Sd = Eb — Ed = 0]. No spin texture forms in the absence of excitonic terms. 
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